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Abstract. Geophysical media are usually heterogeneous and contain multiple species 
of scatterers. In this paper a model is presented to calculate effective permittivities and 
polarimetric backscattering coefficients of multispecies-layered media. The same 
physical description is consistently used in the derivation of both permittivities and 
scattering coefficients. The strong permittivity fluctuation theory is extended to account 
for the multiple species of scatterers with a general ellipsoidal shape whose orientations 
are randomly distributed. Under the distorted Born approximation, polarimetric 
scattering coefficients are obtained. These calculations are applicable to the special 
cases of spheroidal and spherical scatterers. The model is used to study effects of 
scatterer shapes and multispecies mixtures on polarimetric signatures of heterogeneous 
media. The multispecies model accounts for moisture content in scattering media such 
as snowpack in an ice sheet. The results indicate a high sensitivity of backscatter to 
moisture with a stronger dependence for drier snow and ice grain size is important to 
the backscatter. For frost-covered saline ice, model results for bare ice are compared 
with measured data at C band and then the frost flower formation is simulated with a 
layer of fanlike ice crystals including brine infiltration over a rough interface. The 
results with the frost cover suggest a significant increase in scattering coefficients and a 
polarimetric signature closer to isotropic characteristics compared to the thin saline ice 
case. 


1. Introduction 

Natural geophysical media are usually inhomoge- 
neous and contain various species of scatterers with 
different sizes, shapes, and permittivities. Much of 
the early work on multispecies media focused on 
calculations of effective permittivities. A mul- 
tiphase dielectric mixture theory was presented by 
Tinga et ai [1973] with a summary of most impor- 
tant mixture relations derived by many authors. In 
a study of heterogeneous earth media such as 
brine-saturated rocks, it was found that the particle 
shape plays a major role in the effective electrical 
properties characterized by mixing formulas [Wait, 
1983]. Expressions for effective permittivities of 
media with multiphase ellipsoids have been derived 
by Sihvola and Kong [1988], whose results lead to 
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generalized Lorentz-Lorenz, Polder-van Santen, 
and coherent potential quasi-crystalline approxima- 
tion formulas with applications to snow and sea ice. 

The above formulations account for absorption 
loss in the inhomogeneous media. To include scat- 
tering effects, the strong permittivity fluctuation 
theory is used to calculate effective permittivities of 
multispecies media such as dry and wet snow and 
vegetation [Tsang and Kong , 1981a, b; Tsang et aL , 
1982]. For sea ice containing brine and air inclu- 
sions, effective permittivities are computed with the 
strong fluctuation theory at microwave frequencies 
[Stogryn, 1987]. For multispecies media with 
densely distributed particles, effective propagation 
constants are obtained for spherical scatterers un- 
der quasi-crystalline and quasi-crystalline with co- 
herent potential approximations [Ding and Tsang , 
1989]. 

Effective permittivities describe wave propaga- 
tion and attenuation in inhomogeneous media. 
However, to interpret physical characteristics of 
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the media from remote sensing radar data, electro- 
magnetic scattering from multiple species of scat- 
tered needs to be derived. Conventional backscat- 
tering coefficients are calculated with a model of sea 
ice containing spherical scatterers in isotropic me- 
dia with rough surfaces [Tjuatja et al. y 1992]. Pola- 
rimetric signatures of sea ice are derived with a 
model for ellipsoidal inclusions embedded in aniso- 
tropic layered media with rough interfaces [Nghiem 
et al. 9 1993b]. These models consider a single 
species of scatterers (two-phase mixture) and are 
applicable when the inhomogeneities are domi- 
nantly consisted of scatterers of the same types. It 
is shown with a theoretical model of forest back- 
scatter [Lang et al., 1993] that different species with 
different permittivities and structures can cause 
large variations in backscattering coefficients of 
pine trees for vertical polarization at C band. Pola- 
rimetric signatures of a layer of multispecies sphe- 
roidal scatterers are also studied based on vector 
radiative transfer theory [Tsang et al., 1990]. 

In this paper, both effective permittivities and 
polarimetric signatures are derived from the same 
physical description of inhomogeneous layered me- 
dia. The difference with the above-referenced pa- 
pers is that this present work presents a consistent 
model for wave propagation, attenuation, and pola- 
rimetric scattering under the framework of the 
analytic wave theory, which simultaneously consid- 
ers multiple species, ellipsoidal shape, large permit- 
tivity differences, absorption and scattering losses, 
preservation of phase information, multiple interac- 
tions with layer boundaries, and polarimetric signa- 
tures of the heterogeneous media. Effective permit- 
tivities are derived with random orientations of the 
multispecies nonspherical scatterers. The distorted 
Born approximation is applied to calculate polari- 
metric scattering coefficients in the covariance ma- 
trix [Nghiem et al . , 1990], These calculations in- 
volve numerical integrations over Eulerian angles; 
however, analytical expressions are also obtained 
under the low-frequency approximation. Polarimet- 
ric signatures of media with multispecies mixing 
and various scatterer shapes are presented in exam- 
ples for vegetation, snow, ice sheet, and sea ice. 

2. Effective Permittivity 

The effective permittivity of a heterogeneous 
medium is composed of a quasi-static part and a 
scattering effect part corresponding to the first and 


second terms, respectively, in the following expres- 
sion [Tsang and Kong , 1981a] 

eeff=e tf + e 0 [?-i dr -<S>]- , -f dr (1) 

where 7 is the unit dyad, e 0 is the permittivity of free 
space, l g is an auxiliary permittivity, and S is a 
dyadic coefficient. For a medium containing N 
species, effective dyadic scatterer £ eff is given by 


N dir Ctt fin x,y,z 

[£eff]ym = 2 I d y\ #1 da p( a ' P' y ) 

/= I jo Jo Jo 


(0) 
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under the condition |[£ eff ]| jm ^ F In (2), subscript / 
stands for species i, k$ is the wave number in free 
space, r/$ /m is the variance of the scattering fluc- 
tuations (determined by (6) below), and p(a, /3, y) is 
the probability density function of orientations 
given by Eulerian angles a, /3, and y. In this paper 
the random orientation distribution of the multispe- 
cies scatterers is considered. This requires the 
isotropic form of the Green’s function G g [Nghiem 
et al 1990], and the auxiliary permittivity tensor 
becomes e g = e~I. 

The size and shape of scatterer species i are 
described with the Fourier transform of the 

normalized local correlation function. For ellip- 
soids, <J>/|(F) takes on the form 


= 


,p. , 

x My mz 


tt 2 (i + k’ 2 ef x , + k;Mt y + k [%) 2 
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in which £ ix > , and £ iz > are three different 

correlation lengths in the local coordinates corre- 
sponding to the minor, the meridian, and the major 
axes of an ellipsoidal scatterer of species i. The 
local principal coordinate system (x\ y f , z’) of an 
ellipsoid is related to the global coordinate system 
(x, y, z) by Eulerian angles a, /3, and y as shown in 
Appendix A. 

To facilitate the derivation, the following princi- 
pal variances are defined in the local coordinates as 


^i£j’k‘ ~ Sitk’j' — (bf)(£ik' ~ €bk') ( 4 ) 

where /, k’ = x\ y\ z' , subscript si denotes 
scatterer (indicated by .v) species /, subscript b is for 
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the background medium, and the following defini- 
tions have been used 


E »S ix'i^si Eg) 


independent parts in the effective scatterer terms as 
required by the secular elimination condition. Ap- 
plying the symmetry of the ellipsoidal function 
leads to the following results 

Siz'fe. V* tv) 


^ -vi c g 

£() + S jy'{£ S i — £ g) 


£q E Eg) 


replaced by b and species permittivity e si by back- 
ground permittivity e h . In the global coordinates the 
Eulerian transformation rotates the principal vari- 
ances Sjg'fr tO 

I Ifjjklm = & a xjk a xlm E & igx'y' a xjk a ylm E & i£x' z' & xjk a zlm 
E &i£y‘x' ttyjk^xlm E 5 i£y‘y'&yjk&ylm E 5 i£y'z' ^yjk&zim 
E & i£z'x' ^z/kElxlm E <5 i£z’y' ^ zjk&ylm E &i£z'z’ ^zjk^zlm (^) 

where a ^{ s have been defined in Appendix B. 

The condition of secular elimination is used to 
derive e g and S [Tsang et a/., 1981a; Yueh et al 
1990; Nghiem et aL , 1993a, b]. This condition 
imposes £tr)) — 0 where the ensemble average is 
obtained by 

<f(r)>= f dy f dp f da p{a, p, y)T ] 

Jo Jo Jo 


fb£b + ^fsih 'T 

i=i 


In (7) the probability density function of orienta- 
tions for randomly oriented scatterers is 

p(a y p, y) = sin /3 / ( 8 7r 2 ) (8) 

which renders e g isotropic. As shown in Appendix 
B, this orientation distribution together with the 
secular elimination condition cast (7) into 
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(5b) 

C2tt 

e 0 (l + ap i— 


p- 



d<t> 

f= [\ a i 

- tan 1 

v«.] 


(10a) 


Jo 

27re 9 a,\a, 





(5c) 

Siy 1 

fV 

II 

Sh-Uts 


V) 

(10b) 

subscript / 

Si* 


- SiAt* 

, V- 


(10c) 


where principal branch cuts and principal Riemann 
sheets have been chosen for the square root and the 
inverse tangent functions. In (10) the integrations 
over <t> can be carried out numerically and a x and y x 
are defined as 


a t = yf - 1 

cos 2 <t> sin 2 <p\ _1/2 

2 + 2 ) 

7ix' T/v' / 


Then, the average dyadic coefficient ( S ) in (1) is 
determined by 


(S) = ^ f 5i f dy l dp \ dap(a,p,y)f 

i= 1 Jo Jo Jo 


'Six’ 0 0] f 

0 Sp 0 E5^)/ 

0 0 S,v 3 


E^J (5,v E5 fV E5, v )/=5? (12) 

i= i 

where 5^ = £/„•$,•/£/„■. Equations (9) and (12) are 
used to solve for and S’s with an iteration 
method. 

To complete the calculation of the effective per- 
mittivity, let the quantity in the brackets in (2) be 


where ^s have been defined by (5). Elements of S; 
are calculated with the removal of the frequency 


V = *o 2 dk’ [G g (k')] k i®if(k f ) E [$,]*/ (13) 
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In the local coordinates, r/,(r') is a diagonal tensor 
whose elements are 


Vil'f = iv = l if + S if ’ f = X ’> (14) 

where involves the integration of dyadic Green's 
function and correlation function. The result for I iz > 
is 


f 2 ” 

~ I d(j> (^i<p 

JO 77 8 9 


(15a) 
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(15c) 

(15d) 
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yf + £ 
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The symmetry of the ellipsoid is used again to find 
7 /x and /,y , which are simply 

/ivWijrs = M**. €*■. Cv) (16) 


/yW/is V* = iiz'Wix'. Cy) (17) 

To rotate rj,(r') into the global coordinates, the 
Eulerian transformation is applied. Elements of 
rn(r) in the global coordinates are 


V ikl Vix'&xkl T Viy'Gykl Viz'^zkl V ilk (18) 

Variance T in ( 6 ) and quantity 77 in (18) are com- 
bined to define 


x,y,z x,y,z 

£ijm ~~ ^i£jklm V ikl ~~ T! ^ i^u'v’ V uxj (l vxm a wxx 

kj u,x\ w 

&uxj ^irym^wxy ^uxj^vzm^wxz ^ ^uyj^vxm^wyx 
+ ^uyj^irym^wyy T U U yjtl vZm U W y Z + a uz ja vxm a wzx 


After the integration over the Eulerian angles for 
the terms in (19), the effective scatterer tensor for 
species i becomes 


l/eff = J dy j dp j da p(a , , y)lf ] 

Jo Jo Jo 

— 3 igx'x' V ix’ + &i£y'y tr Uy' T $i£z'z 'Viz' V (20) 

Another method to find (20) is to carry out all tensor 
operations involving T’s and rj's in the local coor- 
dinates and then rotate the result into the global 
coordinates; Appendix B is readily applicable to 
arrive at (20). Substituting (20) in (2) and using the 
result in ( 1 ) yield the effective permittivity 


^ eff 


E n 


eo 


X Vix' T ^igy'y' T)iy' + $ii;z'z' V iz‘) 


(21) 

The effective permittivity, as indicated in (21), is 
isotropic due to the random orientation of the 
scatterer species. The above result for the effective 
permittivity of a heterogeneous medium containing 
multiple species of scatterers will reduce to that of 
a two-phase mixing medium containing a single 
species of randomly oriented spheroids embedded 
in a background medium [Yueh et al . , 1990; Nghiem 
et aL, 1993a]. In the next section the effective 
permittivity is used in the derivation of a complete 
set of polarimetric backscattering coefficients under 
the distorted Born approximation. 


3. Scattering Coefficients 

For applications to polarimetric remote sensing 
of geophysical media, it is necessary to derive the 
full set of polarimetric backscattering coefficients 
constituting the polarimetric covariance matrix. An 
electromagnetic wave incident on heterogeneous 
media is scattered by embedded inhomogeneities 
which also dictates the wave propagation and atten- 
uation characterized by the effective permittivity 
obtained in the last section. Thus the same physical 
description of the media used to calculate the effec- 
tive permittivity must be used for a consistent 
derivation of the scattering coefficients. 

Consider an electromagnetic wave impinging 
upon layered inhomogeneous media in the direction 
of wave vector k {)i with an incident angle of %, as 


T ^uzj^vym^wzy T ^ uzj^vzm^ wzz ) 


(19) 
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2 = 0 


Region 1 Heterogeneous Medium 
Ni = number of scatter species 
Itix’iUiyultix* = corr. lengths 
= permittivity of species t 
f $t \ = volume of species t (%) 
c bl = background permittivity 
f b i = background volume (%) 

— 2 = —d\ 

Region 2 Heterogeneous Medium 
N 2 = number of scatter species 
Iflx'Jfiy'Jtw = corr. lengths 
€ 3t 2 = permittivity of species t 
fst2 = volume of species t (%) 

4 2 = background permittivity 

f b 2 = background volume (%) 

z = — 

Region 3 Homogeneous Medium 
c 3 = permittivity of medium 


Figure 1 . Inhomogeneous media with multiple species 
of scatterers. 


shown in Figure 1. The upper half space is air and 
the underlying medium is homogeneous. Each re- 
gion in the layered configuration is inhomogeneous 
and contains multiple species of scatterers which 
are ellipsoidal and randomly oriented. Polarimetric 
scattering coefficients of the layered media are 
defined with ensemble averages of scattered fields 
with all polarization combinations. These averages 
are calculated with spatial integrations over prod- 
ucts of dyadic Green's function G (DGFs), mean 
field F, and correlation function C as follows 


2 Nn 

(E 0s (r) ■ E* 0s (r)) = 22 

*=l t = 1 


2 *o 4 P^yf>P 

d0 JO Jo 


da 



dr*n Pn(a , p. 
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a, P, y) 


• l(Gonij(r, r n ))(F nk (r n ))] • [<G 0 „ l7 (f, 0)(F nm (0)]* (22) 

where t stands for the scatterer species /, and N n is 
the number of species in region n = 1,2. The DGFs 


and the mean fields are given by Nghiern et al. 
[1990]. 

Physical properties of the scatterers are charac- 
terized by correlation function C defined in the 
spatial and spectral domains as 

G igryklmif n i ^ ) 


= <Grvk(r n )tfnim(0\(x(r n ), P{r n ), y(f n )) (23a) 

CtWUmiTn* Ot y p,y) 

= P dp n= 1,2 (23b) 

which are expressed in the global coordinate system 
( x , y, z) and related by the rotation transformation 
with Eulerian angles (a, /3, y) to those in the local 
coordinate system (Jc\ y\ z f ). Function is a 

product of variance T tn j k i m and normalized func- 
tional part $> { £ n 


® tnjkim (P) ~ P tnjkim $ t£n (P ) (24) 

where <t> t £ n is defined in the same manner as in (3) for 
ellipsoidal scatterers with ( ix >, and Z iz > replaced 
by £ tnx «, and t tnz! for species t in region n. 

Variance Y tn j k i m is related to local principal vari- 
ance 8 tr y’k' which is defined as in (4) but the last 
term in the product is complex conjugated 

Stnj'k' = &lnk’f = f stnf bn tnj' ~ hty'Mtnk' ~ Sbnk>)* (25) 

Note that all quantities in the right-hand side of (25) is 
the same as in the calculation of the effective permit- 
tivity; here subscript t (instead of / as in (4)) is used to 
denote the species and subscript n is for region n as 
indicated earlier. For ^s in (25), (5) is used with 
subscript i replaced by tn y e in for the scatterers, e bn for 
the background, and e gn for the auxiliary permittivity 
in region n . After the Eulerian transformations into 
the global coordinates, element Y tn j klm in the fourth- 
rank variance tensor becomes 


*.y,z 

P tnjkim ~ & tnuv&ujk a vlm “ & tnx’x’ a xjk&xlm 

u,v 

F & tnx'y' a xjk a ylm F ^ tnx'z' &xjk a zlm ^ & tny’x' a yjk a xlm 
F $ frty 'y' Qyjk Qylm ^ tny'z' ^yjk&zlm F &tnz’x‘ ^zjk^xlm 

F & tnz'y' a zjk a ylm + $ tnz' z' a zjk a zlm 


(26) 
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When integrations over Eulerian angles of the prob- 
ability density function (8) for random orientations 
are applied to the variances in (26), the averaged 
variance is defined, for latter consideration, as 


in which - 1 or d is for downgoing waves and 1 or u 
for upgoing waves. Expressions for 'P’s have been 
derived by Nghiem et aL [1990], Quantity W^f m 
for region 1 in (29) is 


A 


tnjklm 



w: 


abed _ 
t\jklm 



da p(a, p, y)T njklm I^ d 


(30) 


'P n {a y 0, y)S tnuv a uJk a vim (27) 

Each of the six subscripts u, j\ k , v, /, or m has three 
components x> y, and z; thus there are 3 6 or 729 
terms to be integrated over the Eulerian angles in 
(27), and among those, 540 integrations turn out to 
be zero. The nonzero results are presented in Ap- 
pendix C. 

With the above correlation functions, correla- 
tions of the scattered fields can now be determined 
by (22) to obtain polarimetric scattering coefficients 
defined by [Nghiem et aL, 1990] 


where // \ bcd is the integral of the correlation func- 
tion and the phase in region 1 over the z component 
of the spectral domain and thickness d \ ; the integral 
is carried out under the low-frequency condition as 

/;f rf = dp z p z ) J° dz, e 

/*Q [ g “ it* ab ~ Kcd)dl J 

dz" e np = 2ni<t> lH (0) 

i d , [ Kab-K<d 

(31) 


4 t rr 2 

* TfJLTVK ~~ lim . p, pi (28) 

r-»>oc « KO 

A — ** 


where subscript o denotes incident fields and s for 
scattered fields, subscripts /a, v, t, or k can be 
horizontal ( h ) or vertical (v) polarization, r is the 
range from the antenna to the target, and A is the 
illuminated area. Conventional backscattering coef- 
ficients are a hh = a hhhh , a vv = <r mmn and a hv = 
^hvhv' Backscatter ratios are defined as y = o’ vv la hh 
and e = o- h Ja hh . Correlation coefficient p between 
the horizontal and vertical returns is p = VhhviJ 

(fhhhh 0*vmn.> ) 

All integrations in (22) are carried out analytically 
except the triple integration over the Eulerian an- 
gles which may be intensive in computer time for 
numerical solution. Under the low-frequency ap- 
proximation, an analytical solution, which is com- 
putationally efficient, is obtained and the result is 
expressed as 


With (26), (27), and (30) Wf\j£f m is arranged in the 
following simple form: 

= A , Wm lT d (32) 

Similarly, quantity Wg ? kim for region 2 in (29) is 


^ dy [’ dft P” da p(a,p, yW l 2 Jklm I% rs 
Jo Jo Jo 

(33) 

where the integral of If^ rs is carried out under the 
low-frequency approximation for thickness ( d 2 - 
d\) of region 2 as 



dPr d>,2i(2£p,, P z ) 



dz 2 



Ni - 1,1 x,y,z 

£ 2 S KrJ^t,,n,KlL 

X- 1 a.b.c.d 
Nz d,u JT,V,£ 

+ ni ^2^Tjk^2t *Kjm WftjUm 

/=! p,q,r,s j,k,l,m 


' ii H pq ^ rj )(^2 p U.K py K f y)d \ 


= 27Tf<t>,2f (0) 




(34) 


( 29 ) 


From the results in (26), (27), and (33), it is ob- 
served that is simply 


Wffim = A , 2jklm lT 


(35) 
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z = 0 m 


Inhomogeneous Medium : N\ = 1 

e bl = (1.00 + 1 0.000) eo 

€,n = (3.15 + 1 0.002) €q 
tux' lily 1 Ziw ar€ variable 
Fractional vol /,u =20% 


Underlying Medium: Homogeneous 
€ 3 = (6.00 + i 0.600) € 0 


Figure 2. Parameters to study effects of scatterer 
shapes. 


Note that the above results for the polarimetric 
scattering coefficients are expressed in the scattered 
basis. Both scattering regions in this configuration 
contain multiple species of ellipsoidal shapes, and 
the scattering media are effectively isotropic due to 
the random orientation. For a simpler configuration 
such as a half space or a scattering layer over a half 
space, the results are applicable by just setting the 
appropriate thicknesses. In the next section the 
model is applied to study the polarimetric scattering 
properties from heterogeneous media. 

4. Results 

4.1. Effects of Scatterer Shapes 

In natural media the shape of scatterers such as 
ice grains in snow is particularly subjected to tem- 
perature, medium metamorphoses, or environmen- 
tal variations. The model in this paper accounts for 
various shapes which can be spherical, prolate or 
oblate spheroidal, or ellipsoidal in general. For an 
inhomogeneous medium with a fixed number of 
scatterers and fractional volume, the shape of the 
scatterers is a significant factor on the effective 
permittivity and the scattering coefficients of the 
medium. To study effects of the scatterer shape, 
consider a layered medium such as that described in 
Figure 2. The scattering region consists of an air 
background and embedded scatterers of 20% frac- 
tional volume and permittivity (3.15 T /0.002)e 0 for 
ice at the C band frequency of / = 5 GHz. The 
underlying medium has a permittivity of (6.0 + 
z0.6)e 0 > which is in the permittivity range for soil. 


Table 1. Correlation Lengths and Effective Permittivities 


Shape 

mm 

Gy. 

mm 

Gy. 

mm 

£ effl^ e 0 

Spherical 

0.15 

0.15 

0.15 

1.29 + i 1.95 x 10 “ 4 

Prolate 

0.05 

0.05 

1.35 

1.31 + (2.27 x 10~ 4 

Oblate 

0.58 

0.58 

0.01 

1.33 + (2.63 x 10“ 4 

Ellipsoidal 

1.35 

0.25 

0.01 

1.32 + (2.64 x I0“ 4 


Four different shapes of scatterers are investi- 
gated: spherical, prolate spheroidal, oblate spheroi- 
dal, and ellipsoidal. The fractional volume and the 
number of scatterers are kept the same in all cases, 
and the nonspherical scatterers are randomly ori- 
ented. For the ellipsoidal shape the longest correla- 
tion length is taken to be equal to that of the prolate 
spheroidal shape, the shortest correlation length is 
the same as that of the oblate spheroid, and the 
meridian correlation is chosen by setting the surface 
area (defined by the two longer correlation lengths) 
identical to that of the oblate spheroid. Correlation 
lengths and calculated effective permittivities cor- 
responding to the various shapes are reported in 
Table 1. It is seen that the effective permittivity is 
increasing when the scatterer shape varies from 
spherical, to prolate spheroidal, and then to oblate 
spheroidal forms. Permittivity results for the ellip- 
soids and oblate spheroids are similar due to the 
similar cross sections of the scatterers in the two 
cases. 

Differences in the scatterer shape are also mani- 
fested in scattering coefficients. Copolarized back- 
scattering coefficients for different scatterer shapes 
are shown in Figure 3. Oscillations in the curves of 
the scattering coefficients are due to multiple inter- 
actions with the layer boundaries. The copolarized 
backscattering increases as the scatterer cross sec- 
tion increases from the spherical to prolate, and 
then to oblate spheroidal form. The copolarization 
results for the ellipsoidal and the oblate spheroidal 
forms are similar. For cross-polarized returns, the 
spherical case gives zero value under the first-order 
distorted Bom approximation. Cross-polarized 
backscattering coefficients for nonspherical scatter- 
ers are plotted as a function of incident angles in 
Figure 4. The oblate spheroids provide the highest 
returns while the ellipsoids and the prolate sphe- 
roids have a similar cross-polarization effect due to 
the same size in their longest correlation lengths. As 
illustrated, the ellipsoids behave as a hybrid be- 
tween the oblate and the prolate spheroids. 
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Figure 3. Copolarized backscattering coefficients H — Figure 5. Magnitudes of correlation coefficient p for 
<r hh and V = <r w for different scatterer shapes. different scatterer shapes. 


A common feature in all of the cases for the 
different shapes is the effective isotropy of the 
inhomogeneous medium. For the case of nonspher- 
ical scatterers the medium is rendered isotropic by 
random scatterer orientations with no preferred 
direction. Wave speed and attenuation in the inho- 
mogeneous isotropic medium is the same for hori- 
zontal and vertical waves. This results in correla- 
tion coefficient p between h and v waves with a 
small phase. Magnitudes of p plotted in Figure 5 
show a similar variation as a function of incident 
angles for the different shapes of the scatterers with 
a stronger decorrelation effect from nonspherical 
scatterers. 

4.2. Effects of Multiple Species 

The multispecies model in this paper allows the 
coexistence of disclike and needlelike scatterers. 



0 10 20 30 40 50 


Incident Angles (degrees) 

Figure 4. Cross-polarized backscattering coefficients 
for different scatterer shapes. 


Effects of multiple species are studied with the 
configuration described in Figure 6, which corre- 
sponds to a soybean canopy [Nghiem et al . , 1993a, 
case 9-03-86] where the canopy was modeled with 
only disclike scatterers. The soil surface is rough 
and all input parameters are given by Nghiem et al. 
[1993a]. 

Now, needlelike scatterers are introduced into 
the canopy by allocating part of the total scatterer 
fractional volume to the new scatterer species 
whose permittivity is taken to be the same as that of 
the disclike scatterers. The needlelike form is de- 
picted by a prolate spheroidal shape with correla- 
tion lengths of 0.3 mm and 5.0 mm. Note that the 



Heterogeneous Medium : N\ — 2 

*M = (1.00 -(- i 0.000) c 0 
*•11 = (28.8 i 9.140) Co 
fiir'.y' — 3.2, t\i 2 > — 0.3 mm 
*•21 = (28.8 + i 9.140) c 0 
^21 xV = 0.3, l 2 u' = 5.0 mm 
Volume /,n+/, 2 i =0.88% 

Rough <7, = 1.3, t» = 12 cm 

* = -0.745 m 

Underlying Medium : Homogeneous 

c 3 = (8.76 + i 1.490) c 0 

Figure 6. Parameters to study effects of multiple species. 
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Figure 7. Backscattering coefficients for different mix- Figure 8. Magnitudes of correlation coefficient p for 
tures of multiple species. For each mixture the top curve different mixtures of multiple species, 
is for (Tfjfj , the middle for and the bottom for cr hv . 


volume of an individual needlelike scatterer is much 
smaller than that of the disclike scatterer for the 
given parameters. Three cases are considered here: 
Case 1 has 100% oblate spheroids, case 2 has 85% 
oblate and 15% prolate, and case 3 has 70% oblate 
and 30% prolate spheroids. 

Backscattering coefficients decrease, as seen in 
Figure 7, for increasing needlelike species, which 
has a smaller total cross section. Figure 7 also 
shows that the decrease in the copolarized returns 
is more than that in the cross-polarized return. 
Again due to the effective isotropy, the phase of 
correlation coefficient p is small for the cases under 
consideration. While the magnitude of p for differ- 
ent mixtures is similar as indicated in Figure 8, 
intensity ratio e = c t^Jo ^ is higher for a mixture 
containing more needlelike scatterers, as shown in 
Figure 9. This will result in a copolarized signature 
with a higher pedestal [Nghiem et al . , 1990] indicat- 
ing a stronger depolarization effect. 

4.3. Applications 

This section presents applications of the multi- 
species model to study polarimetric scattering sig- 
natures from inhomogeneous media. A very lossy 
species of scatterers with a high permittivity 
strongly affects wave propagation, attenuation, and 
scattering in a multispecies medium. This model 
includes such a species together with other scatter- 
ing species to investigate effects of moisture in a 
scattering medium such as polar ice sheets. Effects 
of frost flowers on thin saline ice are also studied 
with multispecies mixtures of fanlike ice crystals 
and brine inclusions. 


Temperature changes cause variations in mois- 
ture content of snowpack such as polar ice sheets, 
which are important to global climate because of 
their vast coverage [Thomas et aL> 1985]. It has 
been observed that 07,/, is sensitive to snow mois- 
ture and generally decreases with snow wetness 
[Stiles and Ulaby , 1980]. The configuration in Fig- 
ure 10 is a half space consisting of an air back- 
ground and two species of spherical scatterers: one 
with fractional volume of 0.3 and a relative permit- 
tivity of 3. 15 + /0.0012 for ice grains, the other with 
variable fractional volume and a relative permittiv- 
ity of 65 + /34 for water. Correlation lengths are 
determined from scatterer sizes [Yueh et al ., 1990]. 
Backscattering coefficient cr vv at 5.3 GHz and a 20° 
incident angle, corresponding to the operating pa- 
rameters of ERS 1 SAR (the first European Remote 



Figure 9. Intensity ratio e = cr hv !a hh for different mix- 
tures of multiple species. 






844 


NGHIEM ET AL.: LAYERED MEDIA WITH MULTIPLE SPECIES OF SCATTERERS 



fcot 




y<>. 



Air : to. 

V 


/ = 5.3 GHz 

Heterogeneous Medium : = 2 

€61 

= (1.00 

+ t 
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Ga 

< = .llmm 


Fractional vol f s2 1 = vary 

Figure 10. Parameters of a half-space multispecies 
model for dry and wet snowpack. 

Sensing synthetic aperture radar), is plotted as a 
function of ice grain size for various moisture 
contents from dry to wet conditions in Figure 1 1 . In 
addition to the inverse relation between cr vv and the 
moisture with a good sensitivity, the results show 
that the change in backscatter is larger for drier 
conditions and smaller grains, and ice grain size is 
important to the backscatter. 

Leads in the Arctic ice pack are important to heat 
exchange with the atmosphere and salt flux into the 
upper ocean [Maykut, 1982]. Surface characteristics 
of thin sea ice in leads are dominated by the formation 
of frost flowers with a high salinity [Perovich and 
Richter-Menge , 1994]. This layer is simulated with a 
mixture containing randomly oriented ellipsoids and 
brine infiltration by a volume ratio of three parts ice 



Grain Radius (mm) 

Figure 11. Backscattering coefficient cr m , as a function 
of grain size for various moisture contents. 



<r 61 = (1.00 -I- i 0.000) € 0 

€,11 = (52.5 + i 44.07) 
f l ix' — f i ly' = f i ]*' — .5mm 
Fractional vol f 3 u = 2.0% 

€,21 = (3.19 + i.0014) e 0 

30^2ix , d-5f2iy*,^2iz' = 3mm 
Fractional vol f a2 \ = 6.0% 

Rough a, = 1.0, t s = 10mm 

z = —0.025 m 

Inhomogeneous Medium : N 2 = l 

€62 = (3.19 + t. 0014) Co 

€,12 = (52.5 + t44.07) fo 
Cor length £ 12x f = .236mm 
Cor length t\ 2y * = .660mm 
Cor length l X2z * - .792mm 
Fractional vol j a \ 2 = 10.4% 

z _ — 0.137 m 

Underlying Medium : Homogeneous 

c 3 = (60.7 + x 40.70) c 0 

Figure 12. Parameters of a layered multispecies model 
for frost-covered saline ice. 

and one part brine [Perovich and Richter-Menge , 
1994]. The ice crystals are characterized with corre- 
lation lengths € Ux ‘ =0.1 mm, e Uy ' = 2 mm, an d€ lu > 
= 3 mm to describe the fanlike form. Covered under 
the multispecies medium is an anisotropic layer of 
saline ice composed of ellipsoidal brine inclusions 
preferentially oriented in the vertical direction in an 
ice background [Nghiem et al 1993b]. The inhomo- 
geneities in the ice layer are dominated by brine 
inclusions and this layer is modeled with a single 
scattering species. Model characterization parameters 
are given in Figure 12 for the layered media. Figure 13 
compares polarimetric backscattering results calcu- 
lated from the model with measured data obtained by 
the Jet Propulsion Laboratory polarimetric scatterom- 
eter at C band (center frequency at 5 GHz) during the 
Cold Regions Research and Engineering Laboratory 
Experiment in 1993 (CRRELEX 1993) [Nghiem et al 
1994] for thin bare saline ice. The curves calculated 
from bare the model compare well with the measured 
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Figure 13. Polarimetric scattering results for bare saline 
ice: squares, circles, and crosses are for measured data from 
the 1993 Cold Regions Research and Engineering Labora- 
tory Experiment (CRRELEX); curves are model results. 


data. These curves are for the layered medium with a 
smooth surface and without the covering layer, which 
correspond to the configuration of the bare saline ice 
for the CRRELEX 1993 measurements. 


Then, the covering layer is included in the model to 
simulate a full areal coverage of the frost flower 
formation. Note that the interface between this for- 
mation and the saline ice has some roughness de- 
scribed by a height standard deviation of a s = 1 mm 
and a surface correlation length of = 1 cm. The 
small-scale roughness on the surface is suggested by 
Perovich and Richter-Menge [1994] to serve as nucle- 
ation sites for the frost flowers. Polarimetric scatter- 
ing from the surface roughness is calculated with the 
small perturbation method and included in the total 
polarimetric signatures [Nghiem et al ., 1993b]. The 
curves in Figure 14a indicate that this covering layer 
can increase copolarized backscattering coefficients 
by several decibels while there is a small effect on the 
cross-polarized return. Figures 14b and 14c show that 
the complex correlation coefficient p has a magnitude 
closer to 1 and a smaller phase; this behavior is more 
similar to the isotropic scattering characteristics. This 
is caused by the random orientation of the ellipsoids 
in the covering layer and the nondirectional feature of 
the surface roughness. Figure 15 presents the normal- 
ized polarization signatures a n at incident angle % = 
35° as a function of orientation angle a c and ellipticity 
angle j8 c . [Nghiem et al ., 1990] for CRRELEX-ob- 
served (Figure 15a) and calculated (Figure 15b) bare 
saline ice. The signatures of bare ice from CRRELEX 
data and model calculation compare well. The signa- 
ture of the covered layered medium in Figure 16 (also 
at 35° incident angle) has a maximum in the middle 
(<x c = 90° and / 3 C = 0°) and a lower pedestal 
corresponding to a weaker depolarization effect 
[Nghiem et al 1990]. Parameters in covariance 
matrices corresponding to the above polarization 
signatures are listed in Table 2, which shows a 
good match of all covariance parameters for ex- 
perimental and theoretical results of bare saline 
ice. 

5. Summary 

This paper presents a model for layered inhomo- 
geneous media containing multiple species of scat- 
tered to calculate polarimetric backscattering sig- 
natures. The species are allowed to have different 
size, shape, and permittivity. The strong permittiv- 
ity fluctuation theory is extended to account for the 
multispecies in the derivation of the effective per- 
mittivity with a random distribution of orientations 
characterized with Eulerian rotation angles. 

Polarimetric backscattering coefficients are ob- 
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(a) Comparison of Calculated Backscatter 



Incident Angle (degree) 


(b) Magnitude of Correlation Coefficient p 



(c) Phase of Correlation Coefficient p 



Figure 14. Polarimetric scattering results for bare and 
frost-covered saline ice: curves without the symbol “S” are 
model results for bare saline ice, and those with the symbol 
“S” are model results for frost-covered saline ice. 


tained consistently with the same physical description 
used in the effective permittivity calculation. Under 
the first-order distorted Bom approximation, results 
for the scattering coefficients contain triple integra- 




Figure 15. Normalized copolarized polarimetric signa- 
tures at incident angle 0 O , - 35° of (a) bare saline ice from 
CRRELEX 1993 data and (b) bare saline ice from model 
calculation. 


tions over Eulerian angles, which can be carried out 
numerically. Analytical results are also obtained un- 
der the low-frequency approximation involving 729 
terms for the random orientation distribution. The 
model indicates that nonspherical scatterers with 
larger cross sections give higher effective permittivity 
and stronger scattering effects. The random orienta- 
tion renders the medium isotropic for all scatterer 
shapes; however, the shapes of the scatterers are 
important to cross-polarized returns and magnitudes 
of polarimetric correlation coefficient p. For a mixture 
of oblate (disclike) and prolate (needlelike) spheroids 
the oblate species contributes more to the total back- 
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Figure 16. Normalized copolarized polarimetric signa- 
tures at incident angle 0 q, = 35° of frost-covered saline ice 
from model simulation. 

scatter, and the prolate species increases the cross- 
polarized intensity ratio e = a hv lcr hh . 

The multispecies model allows the inclusion of 
moisture in a scattering medium. Besides a high 
sensitivity of backscatter to moisture content in 
snowpack as observed by Stiles and Ulaby [1980], 
results for backscattering coefficients with vertical 
polarization at C band (ERS 1 polarization and 
frequency) indicate a stronger dependence to the 
moisture for drier and finer snow and ice grain size 
is important to the backscatter. Effects of a multi- 
species covering layer such as frost-covered saline 
ice are also studied. Theoretical results for bare 
saline compare well with measured data at 5 GHz 
from CRRELEX 1993. Then, the model is used to 
simulate a full areal coverage of frost formation 
with a layer of fanlike crystals including brine 
infiltration. The results reveal that the frost cover 
with a rough interface can significantly increase the 


backscatter from thin saline ice and the polarimetric 
signature becomes closer to the isotropic character- 
istics. The strong effects of frost flower formation 
suggest that an experiment needs to be conducted 
to study in detail their impact on polarization sig- 
natures of young sea ice and to provide good 
coordinated backscatter and ice characterization 
data for a better direction and substantiation of the 
sea ice modeling. 

With the consideration of multiple species the 
model advances toward a better representation of 
geophysical media. However, the model needs fur- 
ther development to cope with the complexity of 
natural geophysical media. For instance, particles 
are coated due to disseminated sulphide mineraliza- 
tion [Wait, 1989]. Ice grains in snow are distributed 
in size, metamorphosed, and clustered [ Colbeck , 
1986, 1987]; anisotropy in snow was also quantified 
[Davis and Dozier , 1989], Columnar sea ice is an 
anisotropic multispecies medium containing air and 
brine inclusions, whose phase distributions are 
thermodynamically dependent [Weeks and Ackley , 
1982] and influenced by environmental effects such 
as temperature variations [Arcane et al., 1986]. 

Appendix A: Eulerian Rotations Between 
Local and Global Coordinates 

To relate the local scatterer coordinate system (*', 
y\ z') to the global coordinate system {x, y , z), 
Eulerian angles a, ft, and y are used as illustrated in 
Figure Al. 

The rotation a(0 < a < 2 tt) about the z! axis, j3(0 < 
P < tt) about the y j axis, and y(0 < y < 27t) about the 
z 2 axis are described with the following relations: 



x' 


cos a 

sin a 

0‘ 


x' 

T • 

1 a 

y ' 

= 

-sin a 

cos a 

0 


y' 


z' 


0 

0 

l 


z f 


Table 2. Parameters at % = 35° in Covariance 
Matrices of Bare and Covered Saline Ice 



Vhh’ 

y. 

e , 

Id 

^P. 

Case 

dB 

dB 

dB 

deg 

Bare ice (CRRELEX 

-26.3 

-0.2 

-10.8 

0.85 

-16.3 

1993) 





-10.2 

Bare ice (Model result) 

-25.2 

-0.2 

-10.9 

0.82 

Covered ice (Simulation) 

-21.4 

+0.9 

-14.3 

0.91 

-2.8 


CRRELEX 1993, Cold Regions Research and Engineering 
Laboratory Experiment in 1993. 
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0 

0 

1 

j 


z 2 
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and the resultant Eulerian transformation which 
rotates the global to the local scatterer coordinates 
is 
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cos y cos p cos a — sin y sin a 
T = cos y cos p sin a 4 sin y cos a 
-cos y sin p 

—sin y cos p cos a — cos y sin a sin p cos a 

-sin y cos p sin a + cos y cos a sin p sin a 

sin y sin p cos p 


Tensor #r) in the global coordinates (*, y , z) is 
obtained with the rotation operation 


£(r) = 1 • |(r') • T = T~ 


o o 

0 o . f 

o o £ Z ' 


£xx ^ XV £xz 

- (yx $yy €yz (B3) 

Jzx tzy tzz_ 

Explicitly, the elements of £(f) in the global coordi- 
nates are 

jj xx = ^-(cos 2 y cos 2 P cos 2 a 4 sin 2 y sin 2 a 

- 2 sin y cos y cos p sin a cos a) 

4 4y(cos 2 y cos 2 p sin 2 a 4 sin 2 y cos 2 a 
4 2 sin y cos y cos p sin a cos a) 

+ f z '(cOS 2 y sin 2 p) = £*' 0 *** + + fz' a :xx 


t -(sin y cos y cos 2 p cos“ a 
4 cos 2 y cos p sin a cos a — sin 2 y cos p sin a cos a 

- sin y cos y sin 2 a) — f v '(sin y cos y cos" p sin~ a: 

- cos 2 y cos p sin a cos a + sin 2 y cos p sin a cos a 

— sin y cos y cos 2 a) — f z (sin y cos y sin 2 /3) 

= fx'<ixxy + fy'Ayxy + < B4b) 

£ xz - f v (cos y sin p cos p cos 2 a 

— sin y sin p sin a cos a) 

4 £ v '(cos y sin /3 cos p sin 2 a 

4- sin y sin p sin a cos a) 4 f ,'(-cos y sin p cos /3) 


= Zx'Uxxy + fv’ «>•*>• + (B4d) 

f vv = ft- (sin 2 y cos 2 p cos 2 a + cos 2 -y sin 2 a 
4 2 sin y cos y cos p sin a cos a) 

4 £ v '(sin 2 y cos 2 p sin 2 a 4 cos 2 y cos 2 a 

- 2 sin y cos y cos p sin a cos a) 4 (sin 2 y sin ^ j8) 

= &a X yy + SyOyyy + Z^zyy < B4e > 

£ yz = — fj>(sin y sin p cos p cos 2 a 
4 cos y sin p sin a cos a) 

- £ v ’(sin y sin p cos p sin 2 a 

- cos y sin p sin a cos a) - f z *(-sin y sin p cos /3) 

= fx-flx,* + fy a w* + f* ,fl *y* (B4f) 

Hzx ~ £x' a xzx "*■ (y’ a yzx 

= fc-flxxz + iy-Oyxz + iz-azxz (B4g) 

£zy = £x' a xzy 1 %y' a yzy £z’ a zzy 

= (x-Oxyz + iy-a„z + Zz'thyz < B4h > 

fjj = f,.(sin 2 /3 cos 2 a) + 4y(sin 2 p sin 2 a) + cos 2 p 

~~ fr '&xzz ^ £y‘ a yzz * ^z'^zzz ( B4 i) 

Integrating over the above elements yields 

CllT f IT C2 tt 

{€xx) = dy\ dp I dap(a,p, y)£ xx 

Jo Jo Jo 


= \(£x' + ty’ + fz') 


(fyy) "* 


f2ir Ctt 

= dy 
Jo Jo 


dpi da p(a , p, y)( yy 


= 5 (ft’ + fv' + €z- 


ITT TT ITT 

(£ zz ) = dy dp da /?(a, j3, y)f 


= ^(fx' + ?v' + fz' 
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f 77 
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<f*v> = <{ yx ) = dy 

dp 

da p(a , p, y)£ yx = 0 
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ClTT 

f 77 

fin 

(Lz) = (tzx) = 1 dy 

dp 

da pi a, p, y)$ zx = 0 

Jo J 

'o j 

'o 



(B5e) 


(fzv) ~ 



da p(a , /3, y)£ zy = 0 


(B5f ) 


For generality, subscript / for the scatterer species i 
and b for the background are not included in the 
notation of £y(j' = x\ y\ z f ) which can take on the 
value of for scatterer species / — 1 , 2, 3, . . . , TV 
or £ h j> for the background. The above results are 
used to obtain (9). When £ is replaced by S, (12) is 
obtained. 


Appendix C: Coefficients in Variances 

Polarimetric scattering coefficients derived in this 
paper involve the terms A tn j klm which are integrals 
over the Eulerian angles of the probability density 
function of the random orientation and the vari- 
ances in the global coordinates. A tnjklm consid- 
ered as a variance averaged over the orientation and 
defined by 


A twklm 



dap n {a< 0, y)8 

Uiuv a ujkQvlm 


There are six subscripts associated with coefficient 
a " s in the above equation. Each subscript can be x, 
y, or z; therefore, there are 3 6 = 729 terms in total. 
After being integrated, 540 terms become zero and 
the remaining nonzero terms are presented in this 
appendix. Coefficient a’s, given in appendix B, 
come from the rotation operation on the tensors. For 
instance, consider the following integration of the 
product of the orientation probability and a xxx a xxx 



Pn(®y Pi y )Q X x x a xxx 


+ 4 sin 4 y cos 2 y cos 2 p sin 2 a cos 2 a 

+ 2 sin 2 y cos 2 y cos 2 /3 sin 2 a cos 2 a 

- 4 sin y cos 3 y cos 3 p sin a cos 3 a 

- 4 sin 3 y cos y cos p sin 3 a cos a] = 1/5 

The results of the integrations show that many 
combinations of a's have one of the values: 1/5, 
2/15, 1/15, or -1/30. Actually, the results can be 
written in compact forms due to the symmetry in 
the combinations of the subscripts. By defining the 
averaging operator L as 

Clrr r 7 r Cl-n 

L = I dy I dp I da p n (a , /3, y) (Cl) 

Jo Jo Jo 


the nonzero terms can simply be expressed as 


UVV&UVV 1/3 

(C2a) 

La nn ,a sww = 2/15 

(C2b) 

L&um’Quww — La Uvw — La uvw a uwv 


— Lt(i uww ci VH , w — 1/15 

(C2c) 

^‘d/vw^svw 1/30 

(C2d) 

where subscripts r, s, u, v , or w can be x, 
Following are some examples: 

y, or z. 



da p(a, p, y)a xxx a xxx 


1/5 


(C3a) 



da p(a , (3, y )a xxx a yyy 


= 2/15 


(C3b) 


dy P dfi P 

Jo Jo Jo 


da p(a, p, y )a xxx a zxx = 1/15 (C3c) 


C2tt C TT flTT 

I dy dp 

Jo Jo Jo 


da p(a,/3, y)a yxy a x xy = -1/30 (C3d) 



• [cos 4 y cos 4 p cos 4 a + sin 4 y sin 4 a 
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